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ABSTRACT 


We  consider  the  hydromagnetie  induction  equation  for  an 
unbounded  incompressible  cosmic  fluid  of  finite  conductivity 
and  attempt  to  find  solutions.  A  partial  integration  is 
found  possible  in  a  linear  velocity  shear  or  when  the  magnetic 
Reynolds  number  is  large  compared  to  unity.  To  complete 
this  integration  or  to  integrate  the  induction  equation  when 
there  is  no  coupling  between  the  components  of  the  magnetic 
field  leads  to  the  "substantial  diffusion  equation’".  .  The 
integration  of  this  equation  is  reduced  to  solving  a 
Fredholm  integral  equation.  Several  relations  of  general 
interest  for  a  Lagrangian  description  of  fluid  flow  are 
obtained. 
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I.  Introduction 

An  analysis  of  dimensional  relations  in  conducting 
fluids  (Elsasser,  1954)  shows  that  Maxwell's  equations  for  a 
cosmic  fluid  can  be  written 


(i) 


Vx B  =  p  (£.+k  v* B)  . 

\y  is  the  fluid  velocity,  and  the  other  symbols  have  their 
usual  meanings.  Eliminating  between  these  equations,  we  obtain 
the  hydromagnetic  induction  equation , 

|^-=  Vx^xfi  -v  Vx\/xBt  (2) 

where 

v  =  (ncrr1  (3) 


is  the  "magnetic  viscosity."  We  will  assume  throughout  that 
the  fluid  is  unbounded,  homogeneous,  and  incompressible.  This 
means  that  v  is  constant  and 

V*  V  =  0.  (4) 

We  will  further  assume  a  stationary  flow,  i.e. 

av=o,  V=V(D.  (5) 

Our  purpose  in  these  pages  will  be  to  obtain  integrals  of  the 
induction  equation  under  these  and  certain  other  conditions. 

The  problem  will  be  treated  as  kinematical,  with  \/  given. 

Unless  otherwise  specified  the  coordinate  system  used  is 
cartesian  and  defined  by  the  unit  base  vectors  C  g,  €  g)  • 

Scattered  throughout  the  work  are  results  of  general  interest 
for  a  Lagrangian  description  of  fluid  flow.  The  most 
important  of  these  are  gathered  in  Appendix  I. 

On  taking  the  divergence  of  (2)  one  obtains,  since  the 
divergence  of  a  curl  is  identically  zero, 

| E‘ V  B  =0  (6) 

If  at  any  time  B  vanishes  everywhere  in  space,  than  at  all  times 

7-B=o 


The  usual  proof  of  (7)  given  in  textbooks  on  electrodynamics 
does  not  include  the  induction  current,  <?  V  xfe.  Equation  (2) 
may  now  be  written  in  the  more  familiar  form 

§S  =\7x\^x£+  v.  V23  ,  (8) 

where  use  has  been  made  of  a  familiar  vector  identity. 

In  seeking  physically  significant  solutions  of  the 
induction  equation,  the  solenoidal  character  of  must  be 
kept  constantly  in  mind.  While  it  is  true  that  any  solution 
of  (2)  is  automatically  divergence-free,  this  is  not  so  for 
solutions  of  (8).  In  fact  taking  the  divergence  of  (8) 
shows  that 

§^-B  =  *mV2VB  •  O) 

Although  its  solutions  are  not  generally  d ivergence-free 
their  divergence,  nevertheless,  obeys  a  diffusion  equation, 
so  that  if  initially  divergence-free,  they  remain  so.  We 
conclude  that  (8)  is  equivalent  to  (2)  provided  (7)  is  true 
at  t  =  0. 
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II.  The  Lagrangian  viewpoint 

When  the  fluid  has  infinite  conductivity,  so  that 
v  =  0,  (1-8)  may  be  written 

3^  =  (£•  7)  V  ,  (1) 

where  d/dt  denotes  the  substantial  derivative,  3/ 3t  -H-  \S»  V  , 
of  hydrodynamics.  In  writing  (1)  we  have  used  the  identity 

Vx(ax  b)  =  a(  V*  b)  -  b(  V •cl).  +  ( t  •  V)cl~  ( &♦  P7)  jb  (2) 


and  imposed  the  condition  that  \7*B  =  0.  This  has  the  form  of 
the  Helmholtz  vorticity  equation  of  hydrodynamics,  which 
equation  was  integrated  by  Cauchy  (Brand,  1947),  using  the 
Lagrangian  formulation  of  hydrodynamics.  The  integral 
was  rediscovered  by  Lundquist  (1951  and  1952)  and  applied  to 
the  hydromagnetic  equation  (1) 

Let  (3) 
designate  the  position  of  a  fluid  particle  in  terms  of  its 
initial  position,  T° ,  and  the  time.  The  Cauchy-Lundquis t 
integral  of  (1)  then  is 

B  =  ( B°-  V°)  r  ,  (4) 

where  y  stands  for  the  operator  with  components  ci/cx? , 
andB°  and  B  refer  to  the  field  measured  at  a  particle  when 
its  position  is  and  ,  respectively. 

The  importance  of  a  Lagrangian  description  of  the 
fluid  flow  is  at  once  apparent.  The  particle  trajectories 
are  solutions  of  the  differential  equation 


(5) 


and  the  initial  position,  B°>  enters  by  way  of  the  constants 
of  integration.  The  solution  of  this  equation  has  been 
discussed  in  a  previous  report  (Skabelund,  1955) .  We  note  that 
if  t  is  held  fixed, 

dxi  o 
dxi  "  5x* 
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so  that 


dx^^dXgdXj  = 


,0,0,0 

dx^^dXgdXj 


in  which 


is  the  Jacobian  or  functional  determinant  of  the  transformation 
r°  ->  b  <r°  ,t).  But  the  volume  is  conserved  in  an  incompressible 
fluid,  hence 

J  =  1  .  (6) 

It  may  be  remarked  that  if  is  the  Jacobian  of  the  inverse 

transformation  F B  °  (  Y' ,  t)  then  necessarily  JJ~1=  1,  and  hence 


J”1= 


l1  =i . 


A  verification  of  (4)  bears  out  that  one  must  deal 
with  differential  operations  cautiously  on  switching  between  an 
Eulerian  or  field  description  and  a  Lagrangian  or  particle 
formul_ation.  For  example.  If  &  is  expressed  in  Eulerian 
terms,  i.e.  if  K  is  a  field  variable,  then 


means  +  (  V  •  V  )  g  ; 


whereas  if  g  is  expressed  in  Lagrangian  form,  where  f  is  a 
particle  trajectory,  then  and  t  are  the  independent 
variables,  and  „ 

,*),*]  meana  . 


dg/dt  as  it  appears  in  (1)  refers  to  the  substantial 

derivative  (8) j  but  in  verifying  (4)  one  must  bear  in  mind 

that  B  is  now  expressed  in  Lagrangian  form,  and  (9)  applies. 

In  component  notation  (4)  reads 

3x 

Bi  =  B°j-5x^  ,  (' 


so  that 


which  is  (1) 


“cTE  “  W 


i  _  a 


5t 


^i  =  B°  5^1 

ax°  J  3x° 


a/at  and  a/ax°  commute. 


Note  that  is  independent  of  t  and  that 
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Equation  (1)  is  equivalent  to  (1-8)  with  v=0  only  if 
its  solutions  are  divergence-free,  but  the  remarks  previously 
made  concerning  the  divergence  of  solutions  of  (1-8)  do  not 
hold  for  solutions  of  (l).  In  fact,  speaking  of  £>  as  a 
solution  of  (1),  we  have 


^(V-B)  +(V-  ?)(V-B)  =  0,  (11) 

rather  than  3  B/at  =  o,  as  (1-8)  would  imply.  (11)  follows 
after  taking  the  divergence  of  (1)  and  using  (2)  as  well  as 
the  identity 

\7.  f  CL  =  V(f.  CL  +  Cf>  \7‘  CL  .  (12) 

It  must  be  established  then  that  solution  (4)  is  solenoidal. 

The  necessity  of  doing  this  has  evidently  been  overlooked  in 
previous  accounts  of  the  Cauchy-Lundquist  integral. 

A  well-known  interpretation  of  (1)  is  that  lines-of- 
force  of  B  are  "embedded”  in  the  fluid.  With  this  geometrical 
picture  in  mind  it  seems  clear  that  no  new  lines  are  generated 
during  the  course  of  the  fluid* s  motion,  so  that  div B  if 
initially  zero  would  remain  zero.  Another  way  of  looking  at 
the  problem  is  to  note  that  (11)  has  the  form  of  a  first-order 
wave  equation,  so  that  div  propogates  like  a  wave  in  an 
inhomogeneous  V-field  (Skabelund,  1955).  Thus  div  B  , 
once  zero,  shoul1'  remain  zero. 

Using  (4;)  the  divergence  of  (4)  is 


+  Bo  a2xi 

3xi  ax°  1  dx^x0. 


(13) 


f 


in  which  the  first  terra  on  the  right  is  the  initial  divergence. 
In  s  emi-vector  form, 

V-B  =  7°  *B°  +  b°  7*  ^ 

The  quantities  V  •  8  ,  which  will  also  be  encountered 

in  the  following  section,  must  now  be  evaluated.  It  is  tempting 


to  say  that  o 

dx. 


3xfc3xj 


3x 


c)  k  =  0. 


This,  however,  is  not  true,  for  x^  is  not  an  independent 

variable  and  3/3x^  does  not  commute  with  3/3x°  - a  fact  to 

figure  prominently  in  the  next  chapter. 

To  achieve  the  desired  proof,  we  begin  by  writing 
the  identities  ~,„o  -.0  - 


.es 

<r 

* 

= 

.  SX2  . 
’  3xp  ' 

-  ^ 
3xi 

3xi 

• 

■  1  we  have 

1  cu 

w 

Ih  0 

3x^ 

3x° 

3x 

— -  + 

3x° 

^3  = 

3x.^ 

< 

3x2 

3x° 

axz 

3x° 

*S 

3x^^ 

,  3X2 

ffs  + 

3Xg 

3x3  _ 

3x1 

Sx2 

3x° 

3xs 

3xi 

3x° 

3x^ 

+  3x3 

3X 

3x3 

3x3 

3x^ 

3x° 

3x2 

TT5  + 
3x1 

3x3 

3x° 

Solving  for  3x^/3x^  gives 


^5  = 


Solving  for  Sx^/dx^  and  3Xg/3x°  and  combining  with  the  above 
result  gives 


-  =  (  7x°)  {  7x°)  . 


Repeating  these  steps  for  the  other  two  sets  allows  us  to  write 


sr  _ 


(Vx°)  x(7x°)  =Vx(XpVx°), 


where  (  &<  ,(3,  ^  )  is  a  cyclic  permutation  on  (1,2,3).  The  last 
form  on  the  right  follows  because  the  curl  of  a  gradient  is 
zero  and,  for  any  scalar  and  vector,  Cj>  and  CL, 

\7  x  cp  d  -  (7^  x  &■■¥  \7  x  CL  ,  (19) 

identically.  Finally,  because  the  divergence  of  a  curl  vanishes, 
we  conclude  that 

— ( o  =  0  .  (20) 

Substitution  of  this  result  into  (15^;  reveals  that 

V-B  =  (7°  •  B°  =  o  ,  (21) 

the  last  equality  on  making  the  initial  field  divergence-free. 

Thie  completes  the  proof  that  (4)  is  an  integral  of  the 
induction  equation  for  infinite  conductivity.  Clearly  (20)) 
is  the  Lagrangian  equivalent  of  7*  /  =0,  the  Euler ian 

statement  of  incompressibility. 
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III,  A  Partial  Integral  of  the  Induction  Equat Ion 
a.  Physical  consideration 


In  the  absence  of  dissipation  we  have  seen  that  (2-4) 
is.  an  integral  of  the  induction  equation.  The  form  of  (2-4) 
implies  that  the  field  at  any  time  is  expressible  in  terms 
of  its  initial  configuration  and  the  deformation  of  the  fluid. 
The  physical  reasons  for  this  analytical  form  are  clear  if 
we  picture  the  lines-of-force  as  being  distorted  when  dragged 
along  with  the  fluid.  If  there  is  dissipation  we  must 
deal  with  the  equation  , 


If  *  (  V- 7)  B  =  (  B-  V)  V  «  V  |72  £ 


if =  (6-  V)  v  +  v  ps  B  , 


which  is  equivalent  to  (1-2)  or  (1-8)  provided  its  solution 
is  divergence-free. 

The  preceding  suggests  that  (1)  might  have  an  integral 
of  the  form 

B  =  (B-  v°)  r,  (£> 

A 

where  &  is  a  function  depending  somehow  on  K  and  t.  The 
presence  of  a  dissipative  term  In  (1)  means  the  lines-of-force 
now  slip  relative  to  the  fluidj  B>  »  then,  is  no  longer 
expressible  solely  in  terms ^of  its  Initial  configuration  and 
the  fluid  deformation,  and  £>,  unlike  B°=  B  (y'°0)nius't  have  the 

A  A 

more  general  form  £(K  ,t).  We  expect  £  to  obey  a  diffusion¬ 
like  equation - one  which  will  reduce  £  to  £(  Y°0)  when  u=0.. 

Such  an  equation  is 

v^2  £  ;  (3) 

A  A 

for  when  v=0,  d  B/dt  =  0,  meaning  that  £  is  constant  as  one 
moves  with  a  fluid  particle  along  a  stream  line,  and  hence 

A  A 

B=  B  (  f'  ) .  That  a  solution  (2  and  5)  exists  in  certain 
cases  will  now  be  verified. 
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III.  A  Partial  Integral  of  the  Induction  Equation 
a.  Physical  consideration 


In  the  absence  of  dissipation  we  have  seen  that  (2-4) 
is.  an  integral  of  the  induction  equation.  The  form  of  (2-4) 
implies  that  the  field  at  any  time  is  expressible  in  terms 
of  its  initial  configuration  and  the  deformation  of  the  fluid. 
The  physical  reasons  for  this  analytical  form  are  clear  if 
we  picture  the  lines-of-force  as  being  distorted  when  dragged 
along  with  the  fluid.  If  there  is  dissipation  we  must 
deal  with  the  equation  , 

|f-H-n/-7)B  =  (B-l7)\/-h-v(72^ 

(i) 

or^=  (fc-P)V  +  v72B, 

which  is  equivalent  to  (1-2)  or  (1-8)  provided  its  solution 
Is  divergence-free. 

The  preceding  suggests  that  (1)  might  have  an  integral 
of  the  form 

B  =  <B-  v°)  r,  (2) 

✓s 

where  &  is  a  function  depending  somehow  on  K  and  t.  The 
presence  of  a  dissipative  term  in  (1)  means  the  lines-of-force 
now  slip  relative  to  the  fluid}'  B  >  then,  is  no  longer 
expressible  solely  in  terms  ^f  its  Initial  configuration  and 
the  fluid  deformation,  and  £,  unlike  B°=  B  (y^O^ust  have  the 
more  general  form  £>(  We  e  xpect  Q  to  obey  a  diffusion- 

like  equation---one  which  will  reduce  ^  to  B<  y°0)  when  v=0.. 
Such  an  equation  is 

»72  £  ;  (3) 

A  A 

for  when  v=0,  d  B/dt  =  0,  meaning  that  B  is  constant  as  one 
moves  with  a  fluid  particle  along  a  stream  line,  and  hence 

A  ^ 

£  =  B  (  b).  That  a  solution  (2  and  5)  exists  in  certain 
cases  will  now  be  verified. 
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b.  The  conditions  for  a  particular  integral 


Assuming  (2)  we  have 


dBj  dB  a  dx„  A 
+  B 


a  dxi 


^  V*-  XJ  .  £ 

3T  =  SI  3x*  T  Dj  SI  3x*  » 

but  the  Bx^/Bx^  are  expressed  in  Lagrangian  terms  so  that 

,  Bx.  Bx.  3v, 

d  i  B  i  _  i 


(4) 


Thus  (4)  becomes 


or  in  vector  form 


3S  ~  5S  SxJ  "  3x° 


dEl  dB3  a*i  *  £■  8vl 

as  -  at  Bj  j;o 


Bx 


j 


Consider  the  quantity  (3*  V)V  »  in  component  form  we 
-have,  using  (2), 


3v,  A  9x .  3v,  3v. 

BJ  -a,  ■  \  -o  s,  ■  -o 


(5) 


(6) 


A 

§£=  tffh  v°)  r  *  <B-  v°)  v.  (6;i 


(V) 


or 


r.1) 


<B  •  t7 )  v  =  (  £•  v°)  v  . 

A  / 

Since  3  is  a  solution  of  (3)  by  assumption,  (6  ;  may  be  written 

=  (  B  •  v )  v  +  v  ( v2  B  *  f°)  r »  (8) 


which  has  the  form  (1)  except  for  the  last  term. 


It  remains  to  examine  \/  3  »  subject  to  (2).  One 


obtains 


o  A- 
3  Bk 

Sx^Bx.  “  5xT3x.  ~o  2  3x .  r”%  o  +  --  T“  -..o' 

J  J  J  j  9xk  j  3xj3xk  Bx  jBx  .BX, 


b2b1 


/N  p 

3Bk  3% 


a3 
3  x. 


(9) 


jj  k 


or 

where 


V2  B  =  (Vs  B  -v°)  v  +  2A  +  C 


and  have  components 

BB,  32x  .  A 

A,  =  -*  — ,  C,  =  B. 
1  Bx,  Bx.Bx^  1 
J  J  k 


-3 

3  x. 


k  Bx^Bx^Bx® 


(91) 
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A 

Finally,  then,  (2)  is  an  integral  of  (1)  with  £  subject 
to  (3)  if  and  only  if 


A  =  C  =  0  •  (ID 

Since  the  partials  with  respect  to  the  x^  and  x°  which  appear 
in  and  0^  do  not  commute,  the  conditions  (11)  will  not,  in 
general,  be  met.  In  fact,  these  conditions  will  be  found  to  be 
very  restrictive.  This  is  not  surprising,  for  according  "to 
(2),  (9;),  and  (11),  the  quantity  which  accounts  for 

the  diffusion  of  JB  ,  will  evolve  in  the  same  way  as  B  itself. 
Obviously,  not  many  types  of  flow  patterns  will  permit  this. 

It  will  be  shown  presently,  however,  that  (11)  can  be 
satisfied  under  certain  conditions.  Accepting  this  fact  for 
the  moment,  let  us  see  what  has  been  accomplished.  The  task 
of  integrating  (1)  has  been  reduced  to  that  of  integrating  the 
auxiliary  equation  (3).  The  net  effect  of  our  efforts  has 
been  to  eliminate  the  term  (  B  *  V  )  V  from  (1).  This  is  a 
considerable  simplification,  for  in  (1)  the  components  of 
are  coupled  through  the  term  (  £  ■  V  )  ^/ •  Th-e  se^  (1) 
of  three  simultaneous  equations  in  all  three  B1  has  been 
reduced  to  the  three  independent  equations  of  (3),  in  which 
the  are  not  coupled.  In  this  respect  (2)  is  a  partial 
Integral  of  the  induction  equation.  More  precisely,  this  is 
so  provided  B ,  thus  determined,  has  zero  divergence. 

From  (2)  we  have 


3B.J,  SB  j  dx.± 

3x ,  3x,  3x° 

i  i  j 


*  B, 


s2 
3  x. 


3x .  3x 


i  j 


(12) 


however,  it  was  established  previously  that,  for  an 
incompressible  flow, 

v 


-  °* 


and  hence. 


v-B  =  v°-  8  . 


(13) 


It  can  be  shown  that  for  the  types  of  fluid  flow  which 
satisfy  our  criterion  (11)  A 

v°  •  B  =  o  ^  7°*  £  ( r°,o)  =  o. 
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and  thus  A 

=  0  if  (7«B)t=0  =  0.  (14) 

This  demonstration  must  await  our  discussion  of  equation  (11)  of 
sec.ty'j  however,  rather  than  interrupt  the  argument  there,  the 
proof  is  given  in  Appendix  TV  B.  A 

As  to  the  physical  significance  of  S  ,  note  from 
(2)  that  for  t  =  0 

B  (  r°.0)  =  fe(r°,0)  ;  (15) 

A 

thus  £(  p  ,0)  is  the  initial  value  of  the  field,  £>  .  The 
equation  (3)  governing  £>  is  a  diffusion  equation  for  an 
observer  moving  with  the  fluid.  It  may  also  be  written 

§f  *  (V-7)  B=  v  Vs  8,  (31) 

which  has  the  form  of  a  wave  equation  with  dissipation 
(cf.  sec.  4  and  Tech.  Report  17  of  the  author).  It  describes 
a  -field  propagating  in  an  inhomogeneous  anisotropic 
V-field  and  at  the  same  time  diffusing.  Its  solution  will 
be  the  subject  of  sec.  4. 

c .  Satisfying  the  conditions  -  Linear  velocity  shear 

Suppose  that  the  streamlines  are  straight}  we  may 
without  loss  of  generality  take  them  to  be  directed  along 
the  %-axis.  It  will  further  be  supposed  that  V  does  not 
vary  in  this  direction,  i.e.,  we  assume 

V=  v(x1,x5)  .  (16) 

The  Lagrangian  trajectories  are  given  by 
dx  dxp  dx 

sre  =  °»  xi  =  v»  xi  =  °*  (17) 

whose  ^solutions  are 

Xi=x°,  Xg=x°  +  vt,  x3=x°  .  (18) 

We  have  seen  in  the  preceding  section  that  a 
partial  integral  of  the  form  (2)  exists  only  when  the  quantities 
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of  (10)  vanish.  If  we  require  the  vanishing  of  these  quantities 
without  placing  any  restrictions  on  the  magnetic  field,  a 
necessary  and  sufficient  condition  is  that 


a2 
8  x. 


3x  ,3x? 

J  k 


=  0. 


(19) 


By  referring  to  (18)  we  see  at  once  that 
32X;l  32x2  32x3 

3x^3x°  3x^.3x°  3x  ^  3x° 


It  remains  only  to  consider  32x0/3x  .3x°  and  32x0/3x,3x° 

d  J  1  2  j  o 

The  first  of  these  can  be  written 


(20) 


,2 

3  x2 

-.2  2 

3  x 

3Xj3x^ 

3xj 

3x j3x^ 

(21) 


Now,  the  3x^/3x^.  may  be  found  at  once  from  (18).  In  more 
general  cases,  solving  the  Lagrangian  trajectories  for  the 
Inverse  functions  x°  =  x°(x^.,t)  Is  no  easy  mattery  still 
the  dXg/dXj  can  easily  be  found  in  terms  of  the  3x,</3x°, 
once  the  latter  are  known--~as  they  will  be  from  the 
trajectories.  This  procedure  will  now  be  discussed  because 
of  its  importance  to  a  general  Lagrangian  formulation, 
although  as  mentioned.  It  is  not  essential  here. 

Return  to  equations  (2-14)  of  sec.  2 j  if  they  are 
viewed  this  time  as  comprising  three  sets  of  algebraic 
equations  for  the  nine  3x°/3Xp  In  terms  of  the  Sx^/Sx®,  then 
we  have,  on  exchanging  the  roles  of  K  and  T°  in  (2-17), 

&°=  (  V°xB)x(V°x  ),  (22) 

3x  p  0 


where  (  <*,(3,^)  is  again  a  cyclic  permutation  on  (1,2,3) 
Using  the  above  result,  or  proceding  directly  from 
(18),  shows  that  all  the  32x4/3x 43x°  are  zero  save  four;: 

32v 

t  3Us)2 


,2 

a  x2 

2  i 
3  v 

„2J  11 
8  x2 

Sx-^Sx® 

*  a(x“)2  * 

axsax° 

aS]C2  = 

a2x2 

a2v 

3x33x.° 

3x^3x° 

-  O^O 

ax1ax3 

(23) 
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The  restrictions  imposed  on  the  flow  by  (19)  are  that 

~,2  ,2  a2 

o  v  _  o  v  _  o  v  _  q 

Ov2  ->0,0  _  ,  Ov2  t 

d(x^)  dx^dx^  3(Xg) 

which  have  a  solution 


(24) 


v  =  ax^  +  bx3  +  c  (25) 

The  equi-velocity  surfaces  are  planes  ax^  +  bx°  +  c  =  const., 
parallel  to  the  streamlines.  Thus  v  must  change  linearly  in 
a  direction  normal  to  the  flow  (the  2,-axis).  We  could  just  as 
well  take  this  direction  to  be  the  1-  or  3-  axis  so  that, 
for  example, 

V  =  ( aXg  +  b)  B2  ,  (26) 

which  represents  a  linear  velocity  shear  in  the  3-direction. 
Many  velocity  fields  will  obey  this  condition  locally,  except 
at  points  where  |v|  has  extrema.  Recall  that  has  not 
been  limited  in  any  way. 


d.  Approximate  theory  for  large  magnetic  Reynolds  numbers 


The  possibility  of  satisfying  the  conditions  (11) 
approxmately  will  now  be  considered.  Suppose  the  larges  of 
the  Bx^/dx®  has  an  order  of  magnitude  cT j  we  write 


3x , 


ox. 


<  r  . 


(27) 


Then,  if  X  is  the  scale  of  the  fluid  motion  and  if  the 
motion  is  reasonably  ’'smooth'’'. 


a2 

8  *i 

3x  .3x? 
2  k 


< 


and 


$ 


a3 

3  xi  ^ 

n  ^  7  > 

8x.Sx.3X°  X 
2-  2  k 


whence 


0  =  |C|rO'B^.  . 


(28) 


(29) 
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We  need  next  to  determine  the  order  of  the  dB^/dx^;-  it  is 
not  B/X,  for  thefe-field  has  not  necessarily  diffused  through 
"a” distance  characterized  by  X.  The  ordinary  diffusion  equation, 
aB/at  =  vV2B,  leads  to  a  decay-time  X2/v,  where  X  is  a 
length  characteristic  of  the  dimensions  of  the  conductor.  This 
may  be  regarded  as  the  time  necessary  for  the  field  to  diffuse 
through  a  distance  X.  One  may  say  that  the  field  diffuses 
with  a  speed  v/X.  The  equation  (3)  is  interpreted  to  mean 
that  the  decay-time  of  the  B  -f ield  for  an  observer  moving 
with  the  fluid  is  X2/v,  and  the  speed  of  diffusion  r elative 
to  this  observer  is  of  the  order 


v  _  v 


(50) 


where  Rm  is  the  magnetic  Reynolds  number  (Elsasser,  1954). 

During  the  time  a  fluid  particle  has  moved  a  distance  comparable 
to  the  scale  of  the  velocity  field,  X,  the  B-field  has  diffused 
through  a  distance 

i  =  Vr.  •  (31) 


We  see  that 


whence  also. 


A 

3B 


'k  .  B  RmB 

*  2 


b2b 


'k 


R  B 
m 


A 


A  =|A\ 


R  B  cT 


m 


(72B*7°)r 


^  ~~Z 


A  comparison  of  tdrms  of  the  right  of  (&*)  gives 

C 


A 


A)  -i 

R2 

m 


(v^-f  )rj  ^  . 


(32) 


(33) 


(34) 
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? 


*  i 


,  ] 


Hence,  if  the  flow  is  reasonably  smooth  both/\  and  Cmay  be 
neglected  in  comparison  with  the  first  term  in  (9*) 
provided 


R. 


m 


£ 


10 


(35) 


The  divergence  condition,  (14),  is  satisfied  to  the  same  order 
of  small  quantities,  as  discussed  in  Appendix  IV  B. 


>  i 

h  f 


K  i 


— r 
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IV.  The  Substantial  Diffusion  Equation 

a.  Prelim inaries ;  nature  of  the  problem 


Suppose  there  is  no  interaction  between  the  components 
of  f^)  ;  this  means  that  ( B  •  V  )  V  =  0,  and  the  induction 
equation  (3-1)  is  reduced  to 


If  +  (V- V)B  =  V  (72b 
or  =  v  72  £  . 


(1) 


This  is  precisely  the  substantial  diffusion  equation  encountered 
in  sec.  3.  It  arises,  in  three  cases:  it  is  the  induction 
equation  when 

(0*  V)  *  =  0  i  (2) 

it  is  approximately  the  induction  equation  for  large-scale 
velocity  fields,  for  when  V  varies  slowly,  that  is  when 

(3) 

so  that  it  may  be  dropped  from  equation  (1-8);  it  is  the. 

A 

equation  for  £  in  the  partial  integral  of  the  complete 
induction  equation  when  neither  (2)  nor  (3)  is  true. 

We  propose  in  this  section  to  obtain  solutions  of 
(1).  The  symbol^  will  be  used  throughout,  with  the 
understanding  it  may  be  replaced  by  q  when  appropriate. 

Whenever  the  solutions  obtained  are  intended  to  represent 
the  complete  magnetic  field  and  not  jus.t  3  »  they  are 
subject  to  either  (2)  or  (3),  and  this  must  be  verified  for 
particular  solutions.  The  field  is  -continually  spreading 
both  by  diffusion  and  by  its  being  (partially)  dragged 
with  the  fluid.  It  is  not  difficult,  however,  to  construct 
special  V-  and  configurations  which  because  of  their 
geometries  satisfy  (2)  at  all  times.  If  (B  *  V  )  V  is 
merely  small,  Initially,  then  (3)  cannot  hold  indefinitely; 
nevertheless,  it  should  remain  true  for  times  of  the  order  of 
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the  period  of  fluid  motion,  v/X*  This  is  because  in  the 
range  of  Rm  where  both  diffusion  and  induction  are  important, 
the  spreading  effects  of  diffusion  and  "differential  dragging" 
tend  to  oppose  one  another.  (This  assertion  i^  based  on  the 
familiar  "'entropy"1  arguments’  that  the  general  effect  of 

diffusion  is  to  smooth  out  inhomogeneities  in  the  field - in 

this  case  inhomogeneities  caused  by  the  fluid  motion.  See 
also  Inglis  (1955)  who  cites  an  example  to  show  the  extent 
to  which  magnetic  field  lines  can  be  stretched  by  "dragging" 
before  their  "slipping"  velocity  relative  to  the  fluid  equals 
the  fluid  velocity.  Lundquist  (1952)  also  gives  an  example. 

It  has  already  been  remarked  that  the  substantial 
diffusion  equation  may  be  regarded  as  a  first-order  wave 
equation  with  dissipation.  A  closer  examination,  however, 
reveals  that  (1)  is  an  equation  of  the  elliptic  rather  than 
hyperbolic  type  (Sommerfeld,  1949,  p.  38  )  ,  and  mathematically 
it  is  more  akin  to  a  diffusion  equation  than  a  wave  equation. 

It  is  in  this  direction,  therefore,  that  we  look  for  a 
solution.  Physically  the  equation  represents  a  combination 
of  propogation  and  diffusion  in  an  inhomogeneous  anisotropic 
"medium";;  consequently,  it  will  not  be  surprising  to  find 
its  solution  cumbersome. 

b .  The  substantial  diffusion  equation  as  an  inhomogeneous 
diffusion  equation 


Trying  to  solve  the  equation,  even  for  simple  fluid 
flows,  dispels  any  notion  of  there  existing  bounded  solutions 
with  v  ariables  separated.  A  Fourier  analysis  and  an 
application  of  the  WKB  approximation  for  slowly  varying  N/ 
were  also  unsuccessful  in  obtaining  physically  meaningful 
solutions . 

We  therefore  choose  to  look  upon  (1)  as  an  inhomogeneous 
diffusion  equation, 

Vs  B  -  I  |f  =  (l1) 

From  this  point  of  view,J5has  "sources"  -v-'5'  (  n/*  V' )  £> 
which  arise  from  its  interaction  with  the  fluid. 
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Using  Green's  function  for  the  diffusion  equation  one 
can  show  that 


£(  r, t>  =  -v"1jv*Jd3r/G(  y-  v')  S' 


(4) 


+v 


«  V 

-1  /i5r'«>  e'> , 

’’V  o  s 

in  which  £  is  the  field  within  a  volume  V  bounded  by  a  surface 
s»  G(  IT  »^|  ^  ^t  )  is  the  appropriate  Green's  function,  ' 
and  denote  y  (  and  A(  ^  t ty  )  respectively,  and 

p 

t+  means  t  +£  ,  f*0.  The  details  of  this  operation  will  be 
found  in  Appendix  II. 

Assuming  that  G  <-v>  and  B  <•«  l/t'  ,  at  least,  as 
the  surface  integral  vanishes  as  S  receeds  to  infinity,  and 
in  the  infinite  domain 

-1  fat  fa3  r'  afv7-  v'  >  s' 


£(  r,t)  =  -i 


-HV 


-1 


a W  £')t,e 


(5) 


It  follows  that 


(  v-  V)h=  -v"1  fdt'  [  d3^  V-7G)  (  v"-  7) 


~1  fd 3r'A  ^ 

J  u 


t*o 


(6) 


The  last  integral  is  a  known  function  or  rather  functional, 
in  terms  of  £(T,0),  which  will  be  designated  by 


T  (  r  =  v 


rl 


(d3r'[v 7g(  r,t|  r',v]8(  r'.°>* 


(V) 
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Then  (6)  has  the  form  of  a  Volterra-Fredholm  equation  for 
(  V  •  V  )g#  viz.  ^  (8) 

(v-  P)8  =  n  -H-  Jdt'Jd3^ '  k(  y,t  jr't'x  v'.  pig', 

0 

whose  kernel  is 

K(f,t  p',t  )  =  -v“V(f)  *  P&(  t  .  (9) 

At  this  point  it  becomes  necessary  to  state  something 
more  about  the  Green’s  function.  The  appropriate  Green’s 
function  for  the  infinite  d omain  can  be  shown  to  b e  the 
impulse  function 

G(R  Jt')  =  vU  M  (4Tivr)"n//2exp(^.R2/4vr)  (10) 

r  =  |  s  T  -  t-t  , 

where  U  is  the  unit  step  function 

I  ■ 

u(r)Jpfo,  t<o 

'"/I,  T>  0. 

(Morse  and  Feshbach,  1953,  p.  894)  Certainly,  theny  our 
assumption  about  the  behavior  of  G  at  ©o  was  justified. 

Furthermore,  since  G=0  for  all  t/>t,  the  time  integration  in 
(8)  can  equally  well  be  extended  over  all  t,  so  we  may  write 

oc  (H) 

oO  r 

(  v*  v  )£  =  n  ^t)  +  Jdt7  jd5r/  k(  r  ,t  | Kt' )  ( i/-  pi  8* , 

a 

which  is  an  inhomogeneous  Fredholm e  quat ion  of  the  second 
kind.  The  kernel  is  asymmetric? 

T  ..  0(T)V.  R  exp(.RS/4vr) 

K - -  (12) 

2v  V  (4nv  <1*  )  ■' 

K  =  r-r'. 

Fredholm  e  quations  for  a  function  of  more  than  one 


* 


it 

% 


h 

k 


S' 


-20- 

variable  are  seldom  discussed  in  the  textbooks j  however,  they 
can  be  solved  using  straightforward  generalizations  of  the 

I 

standard  methods  used  on  equations  for  a  function  of  one 
variable  (Courant  and  Hilbert,  1953,  page  152) •  One  method 
of  solution  is  found  in  Appendix  IV.  Finally,  let 

f  (r.t)  =  (V-F)B  (13) 

be  the  solution  so  obtained;  then  substitution  Into  (5)  gives 


oo 

g(  r,t)  =  Jd3  r,t  | r1  ,*')  $ 

o 


(14) 


+v 


B(r 


It  is  proved  In  Appendix  IVA  that 


\7*£)(^,t)  =  0  provided  P*£(  p, 0)  =  0  . 


(15) 


c .  The  substantial  diffusion  equation  reduced  to  an 
Inhomogeneous  Helmholtz  equation’ 


Examination  shows  that  there  exist  no  physically 
meaningful  solutions  of  (1)  with  variables  separated. 


(f,t)  =  X^(x,1)Ye<(x2)Z;(x3)T(t) 

(no  summation  on  greek  indices).  This,  however,  does  not 
preclude  the  existence  of  solutions  composed  of  a  sum  of 
terms,  each  term  of  which  has  variables  separated.  In  an 
effort,  then,  to  find  solutions  of  (1)  less  awkward  than 
(14)  we  look  for  solutions  of  the  form 


B(  i\t) 


ik  c  (k)T(t  jjk)  p(  p  ;k) , 


t 


(16) 
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a  superposition  of  terms  with  variables  separated. 
Then 


vp2e-  |f  -  (  /♦  7>B 


=  ^dkc(k)  [t(v  v2f-(v  ?)F -F&] > 


which  vanishes  if 


(17) 


f  g*vA-0 

-vk2f 

(18a)  gives  T  =  e  , 

as  us  ual . 


(18) 


(19) 


The  plan  now  is  to  look  upon  (18b)  as  an  inhomogeneous 
Helmholtz  equation.  Doing  this  one  may  show  that 


F ( r }k)  =  -v‘ 


d 


)f' 


*  ~F/  §£');ds/  •  (20) 


The  details  of  this  operation  are  contained  in  Appendix  III. 

G  is  the  Green1 s  function  for  the  Helmholtz  equation  and  may 
be  taken  as 

G(r/f/jk)  =  R"1exp(ikR)  (21) 


(Morse  and  Feshbach,  1953,  p.  891).  As  S  recedes  to 
infinity  the  surface  integral  vanishes  if  F  goes  at  least  as 
r  as  r  +  .  Assume  this  is  so.  In  the  infinite  domain, 

then , 


F  (  r$k)  =  ~v 


d5r/G(rjr^k)(v/.\7/)F/ 


(22) 


This  is  a  particular  solution  of  (18b) j  to  it  we  may  add 
any  solution,  ,  of  the  homogenous  Helmholtz  equation. 
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72p  +  k2f  =  0  . 


Doing  this  will  guarantee  our  obtaining  a  solution  F  with  zero 
divergence  (a  subject  discussed  in  Appendix  IVC),  We  have  now 

v  ’  r 

p(r;k)=F(  i-k)  »v_1  d5r/G(r)r/ik)(y/*7/)F^  •  (24) 


It  follows  at  once  that 


(  V«7)P  =  (V-V)F  -V1Jd;3r/(  V-Vg)(V/'V/)F/ ,  (25) 

which  is  an  inhomogeneous  Fredholm  equation  of  the  second 
kind  for  (\/*\7)F>  with  the  asymmetric  kernel  -v"1  V*  P'g. 

It  can  be  solved  using  extensions  of  standard  methods.  Let 

/\(  r  °>  k)  =  (  P7)  F  (26) 

be  the  solution  so  obtained;  putting  this  into  (24)  gives 


F(  r?k)  =  f(r  jkj-v"1  Jd9  r;G( rjr^k)  A  ( r) 


k)  .  (27) 


The  general  (though  not  necessarily  complete)  solution 
of  (1)  is 


B  (  r»t)  =Jik  c(k)  F (  Y'  ;k)  exp(-vk2t). 
Furthermore, 

B(r,o)  =Jjkc(k)F(r?k), 


which  is  a  Fredholm  equation  of  the  first  kind,  with  the 
kernel  F  (|*;k).  Its  solution  (if  one  exists),  obtained  by- 
standard  methods ,  gives  the  spectrum  function  c(k)  in  terms 
of  the  initial  field  £  (  ^,0).  As  to  the  divergence  of  £  , 
note  that 

V9  B>  c(k)  (  V«  F)©xp(-vk2t)  , 

which  vanishes  identically  only  if  Vmfr=  °«  In  Appendix  IVC 
it  is  demonstrated  that 

0  provided  VmF=  0  .  (j 
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The  order  of  the  definite  integrals  involved  has  been 
reduced  by  one  by  using  this  quasi-separation  of  variables  in 
(16).  It  has  been  done,  however,  at  the  expense  of  assuming 
that  there  exist  solutions  of  (1)  of  the  form  (16)  which  are 
regular  at  ©oand  which  allow  the  spectrum  function  to  be  determined 
through  (29).  Whether  or  not  this  is  so  must  evidently  be 
decided  in  individual  cases.  In  obtaining  solution  (14)  its 

form  was  in  no  wise  restricted;  hence  there  can  be  no  doubt 
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that  it  vanishes  at  least  as  r  at  infinity,  whenever  the 
initial  field  is  confined  to  a  finite  region (because  from  a 
distance,  then,  the  field  appears  to  arise  from  multipole 
sources ) • 

It  seems  now  natural  to  inquire  if  (1)  could  not  be 
treated  as  an  inhomogeneous  wave  equation  to  which  the 
familiar  Kirchhoff  method  of  integration  might  be  applied. 

This  method  is  not  applicable,  it  turns  out,  because  of 
the  intrinsically  anisotropic  character  of  the  propagation 
which  (1)  describes. 

d.  The  solut ions  in  Lagrangian  form 

Consider  the  subject  of  section  4a.  Once  the  solution 
^  of  (11)  Is  known  one  may  proceed  directly  to  find 3  in 
Lagrangian  form.  Using  (2-3),  (2-5),  and  (13)  we  may  write 

(si) 

where  the  use  of  the  total  derivative  means  we  are  thinking 
of  as  expressed  in  terms  of  Lagrangian  variables.  Then, 
since„d^=  dr*  Vcf  is  the  change  in  in  the  direction  of  df*, 

dB  =  (dr •  7)B  =  $ [r(  r°»*)»t3dt  (32) 

gives  the  change  in  &  arising  from  a  displacement  dr  along  a 
Lagrangian  trajectory  (and  includes  the  effects  of  bo th  time- 
and  space-changes) .  As  a  result 

t 

&(  r,t)  =B(r°so)  *  Jf  ( r°»t)dt  . 


(33) 


The  Lagrangian  and  Eulerian  formulations  are  merely  different 
ways  of  expressing  the  same  thing,  and  because  of  this  the 
conclusion  (15)  about  divB  is  still  valid,  and 

V'2>  =  0  if  7°«  8(r0,°)  =  0  (34) 

Although  it  is  usually  more  convenient  to  have 
solutions  expressed  in  Eulerian  terms,  it  may  be  desirable  to 
settle  for  a  Lagrangian  form  and  thus  avoid  the  complicated 
multiple  integrations  involved  in  (14). 

The  solution  which  is  the  subject  of  section  4c 
may  also  be  expressed  in  Lagrangian  form,  at  a  considerable 
saving  of  labor.  From  (26)  we  have 

dF  =(d  fV)F  =/\  (V (  .  dt, 

t 

and  f~(r)  =  f(r°)  +  J/\(r°,t)dt  . 

Besides  requiring  that  =0,  as  in  (30),  it  is 

necessary  to  have 

V°'F'(r°)=0>  (37) 

because  the  term  F(^°)  is  arbitrary  as  far  as  the  Lagrangian 
integration  is  concerned. 


(35) 

(36) 


)ure-dif fusion  solutions 


We  close  this  discussion  with  an  account  of  one 
class  of  solutions  of  (1)  available  at  once.  This  class  i's 
obtained  by  setting 


St  " 


vf2  JB 


(/•?)£=  o  . 


(38) 


B  obeys  the  conventional  diffusion  equation  and  in  addition, 
according  to  (38b),  must  not  vary  in  the  direction  of  V. 

From  (38)  we  see  that  V*  8  likewise  obeys  a  diffusion  equation 
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so  that 


7*8  =  0  if  (  7*  8  )t=0  =  °*  (29> 

Eet 

XT1(r)  =  k±,  TJg(r)  =  k 2  (40) 

be  linearly  independent  (  x  7ug  =^0)  and  integrals  of 

d  p  x  V=  0  ,  (41) 

This  means  the  intersection  of  the  surfaces  defined  by  (40)  falls 
along  the  streamlines  $  and  if  U3  =  K  x  7bg,  we  have 
(IJj,Ug,U3)  forming  a  system  of  curvilinear  coordinates.  As  a 
consequence  (38b)  is  satisfied  if 

2 

8  =B(U1,Ug,t)  =  e"v  t  r(VlfU2).  (42) 

The  solution  of  (38a)  is  now  a  question  of  separability.  The 
space  part  of  the  solution  satisfied  the  vector  Helmholtz  equation, 
which  for  other  than  cartesian  coordinates  should  be  written 

Y  x7 x  p  V' F~  =0.  (43) 

Its  separability  is  essentially  limited  to  six  different 
coordinate  systems  (Morse  and  Feshbach,  1953,  p.  1767);  viz. 
circular,  elliptic,  and  parabolic  cylinder  coordinates  $ 
spherical  and  conical  coordinates;  and,  of  course, 
rectangular  coordinates.  Thus  equation  (38)  can  be  solved 
simultaneously  if  the  streamlines  of  V  coincide  with  one  set 
of  coordinate  lines  for  any  of  these  six  systems.  If  they 
coincide,  say,  with  Uplines  (the  lines  on  which  and 
Ug  are  constant)  then  p  =  p  (U^Ug)  may  be  any  divergence-free 
solution  of  the  two  dimensional  Helmholtz  equation  In  and 
Ug.  The  solution  of  (43)  Is  completely  treated  elsewhere 
(Morse  and  Feshbach,  Vol,  II,  pp.  1762-1767).  The  simplest 
example  that  comes  to  mind  is  where  g 

y=  e1vand8=e"v  t.P(Xg,x3) 
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with  t-8  J  -g  *  k2)F  =o,  V-F  =  o. 

ax2  Sx3 

These  pure-diffusion  and  non-wave-like  solutions  are  of 
limited  physical  interest.  To  require  in  addition  that 
( e>  •  v )  v  =  0  would  be  altogether  too  restrictive y  hence 
they  are  of  interest  not  as  integrals  of  the  induction 

✓A 

equation,  but  as  possible  -fields,  with  e  given  by  (3-2), 


f 


\ 


f 


-•/Wt**:: 
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Resulta  of  General  Interest  for  a  Lagrangian 
Description  of  Fluid  Flow. 


Let  Y 0  locate  the  initial  p osition  of  a  point  moving 
with  the  fluid.  Its  position  after  a  time  t  is  given  by  the 
trajectory  Y  =  J*  (  »  where  r  is  a  solution  of 

dr  _  v 
dl  "  V' 

\/  being  the  fluid  velocity.  In  the  Lagrangian  formulation 
r°  and  t  are  the  independent  variables  and  d/dt  and  d/dt 
commute  with  the  d/dx°j'  however,  the  d/dx°  and  d/dx^  do  not 
commute.  If  operating  on  a  function  expressed  in  Euler ian 
t  erms ,  then  df/dt  means  d  F?dt  +  (\/»\^)p,  the  familiar 
substantial  derivative^  but  if  operating  on  a  function 
expressed  in  Lagrangian  terms,  dp  /dt  means  simply  d  P  /dt. 

The  following  results  hold. 

J§T=  (V°  XTo)  x  (,V°?C)  =  v°^(4v0xy) 


where 


is  the  Jacobian  or  functional  determinant  of  the  transormation 
r°~>  p(T'°,t),  and  (  o(,p ,p)  is  a  cyclic  permutation  on 
(1,2,3)  .  If 


denotes  the  Jacobian  of  the  inverse  transformation 
.  r-*r°(  r,t),  then  similarly 

j-1  2T  =  (  vx2)  x  (  Vx°)  =  V x  (x°.  Vx°) 
3x*  «  o 

When  the  flow  is  incompressible 

V*  V  =  0 
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and  J  =  J'1  =  1, 

(7) 

which  means  that 

7-  =  v°  • 

i  =o. 

(8) 

j 

(8)  is  the  Lagrangian  equivalent 

of  the  Eulerian  statement  (6) 

expressing  incompressibility. 

Equation  (1)  can  be  used 

in  special c  ases  to  derive 

conservation  theorems  for  T*  .  Suppose  0  denotes  any  linear 
differential  operator  with  terms  in  the  8/8x°.  Then,  using 
our  commutation  rules 


&X,°r-£°v, 

=£°  V=  o 

whence  &C°T  -  const.  =%s°V°  for  any  fluid  particle.  By 
way  of  example,  if  in  particular  cases  V0  • 'V  =  0,  it 
follows  that 

V°  •  r  =  7°*  T°  =  3. 

For  such  purposes  it  may  be  convenient  to  express  \7° 
in  any  of  the  equivalent  forms 


and  if 

we  would  have 


.ii 

8Vi 

=  *_r 

ax° 

a  ° 

°xi 

3XJ 

8x^ 

*  V  vi  =  ^  V1  •■(V3Cg)  x  (Vxj) 


etc. 


Appendix  II.  Integration  of  the  inhomogeneous  diffusion  equation 

The  substantial  diffusion  equation  of  sec.  4  is  to  be 
written  as 

v2b  -  v-1 

and  regarded  as  an  inhomogeneous  diffusion  equation.  Let 


(1) 
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(1)  be  multiplied  by  a  function  G  which  is  essentially  an 
integrating  factor  and  remains  to  be  determined.  Then  let  the 
quantities  3 and  v"^  9S/9t  be  both  added  and  subtracted  on 
the  left.  After  rearranging  terms  and  integrating  there  obtains 


1$(  TT  )  and  \/(  T/,t/  )  are  r  epresented  by  and  V  ,  and 
t  means  t  +■£  ,  the  volume  V  remains  as  yet  arbitrary. 

Now,  let  G(  K  ,t  |  T^t7  )  be  the  response  of  the  system 
at  (.  r,t)  to  an  instantaneous  point  source  of  unit  strength 
located  at  (Tf^t*),  i.e., 

•\72G-v“1|^  =  -  cl(T  •  T)  cT(t-t  );.  (3) 

and  require  that 

&(  r, t)  rW )  =  o,  t<t;  (4) 

(this  is  a  causality  statement,  required  by  the  unidirectionality 
in  time  implicit  in  the  diffusion  equation).  It  is  a  consequence 
of  these  t  wo  conditions  that 

IrU' )  =  gc  r',-t'\r,-t),  (5) 

(See  e.g.  Morse  and  Feshbach,  Vol.  I.,  p.  858)  which  is  a 
statement  of  reciprocity  and  is  rather  obvious  from  physical 
considerations.  It  follows  that 

y2G  +  v"1  22  =  -  cT (  T-  Tf)  <P (t-t  ).  (.6) 

at 


Returning  to  (2)  and  using  this  result  gives 


rdtj ds r '  [\712 g  v“1^/jB(  tW ))  =  - 


As  to  the  last  integral  on  the  left  of  (2),  reverse  the 
order  of  integration  and  obtain 

(8) 

-v^jdV '  j'c it;  2-,  (s  )'  =  -v"1  r,t  |rU'  >3(  rU')j 


/  0 

But  according  to  (4) 


g<  r,t  |  r',t+)  =  o 


(which  was  the  reason  for  integrating  to  t+) y  hence, 

(8)  =  v_1  fd3r  £>')  ,  • 

J  t=0 

Finally,  we  apply  the  symmetrical  form  of  Green's  theorem  to 
the  first  integral  in  (2)'j 


^  r 

jdt/  d 5r  [g  F'V-  e!  V/2g] 

J 


o  V 


*  S 

after  which  (2)  becomes 


"T  ^ 

BCr^)  =  "V-1  Jdt;  Jd3rG(  v"1  d3r/(G^/)'t  =c 


•  / 


at'  ([g  —  ,1  ds/ » 

L  3n'  3n; J 


J  JL  an'  ~an'J 
0  £ 

which  is  equation  (4-4). 

(12)  has  a  simple  physical  interpretation.  The 
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contribution  to  3  from  "’sources"  exterior  to  the  region  enclosed  by 
S  is  expressed  by  the  surface  integral.  The  contribution  from 
within  V  comes  from  an  integration  over  a  continuous  distribution 
of  sources  having  strengths  -(\/»  7)B  and  3  (  Y *0) »  the 
latter  are  instantaneous  s ources . 

The  essential  steps  in  this  integration  may  be  found  in 
Morse  and  peshbach(1953,  pp.  859  ff.)  or £1955,  Chap  .  VII}* 

Appendix  III.  Integration  of  the  Inhomogeneous  Helmholtz 
Equation . 


Equation  (4-18b)  reads 

V8F  *kzF  =  v-1!  v-V)F,  (1) 

a  H.iinholtz  equation  with  "sources”  -v"^(  \S  •  ^7 )  P  •  Let 
G  be  the  response  to  a  unit  point  source  at  T*',  then 

*  k2&  =  -  cT (r-  (2) 


G  may  be  used  as  an  integrating  factor.  Multiply  (1)  by  G 
and  (2)  by  P,  subtract  the  two  equations  and  integrate; 
one  obtains 


Ja3r'(  fWs  v1*  F1  )  (3) 

-  -  fd3r/  cT(  r-  t')  F(  rVv^JdVo  (./•/)  F  • 

V  V 

On  using  Green’s  theorem  in  its  symmetrical  form, 


F  <r) 


-V 


-1 


c 


J 


d5r/  g(  v •  v)  F ' 


-H 


r 


(F 


j 

Z 


f  dG 
dn1 


G 


(4) 


This  is  equation  (4-20) 
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Appendix  IV.  On  the  Divergence  of  Solutions  of  the  Substantial 
Diffusion  Equation. 

A.  The  substantial  diffusion  equation,  repeated  from  sec.  4  is 

|f  *  «V  ■  V\B  =  *•  72E  ,  (i) 

in  which  the  symbol  Brnay  represent  also,  when  appropriate. 

Prom  (2-2)  and  (2-12)  we  see  that 

v  .  V)b=  v-  Vi  y-B)  *  V*  (  8  *F  )  \z“  (2) 

Assume  for  the  moment  that 

(fc-V  )  =  0  or  is  negligible,  (3) 

which  must  be  the  case  if  is  an  integral  of  the  induction 
equation,  according  to  section  4a.  On  taking  the  divergence 
of  (1), 

Jg(  )  4f  V-  V(  V’B)  =  vy2(  ^E),  (4) 

so  that  B  and  its  divergence  satisfy  the  same  equation. 

Let  for  a  moment  the  symbol  b  represent  both  and 
dlv  B  .  WJe  then  write 

V2fc-v-1||3=v-1(  V-  V)b  (5) 

and  discuss  both  relations  simultaneously. 

Suppose  the  "’source"'  function  on  the  right  is 
initially  zero,  then  momentarily  (5)  reduces  to  the  homogeneous 
diffusion  equation. 


dfc  _  xr  2  L 


which  has  solutions  of  the  form 


\ 


I  1  l 


\  I 
)  f 


I  >  ! 

•  ■  i .  t 
:  i 

!  J 


i  l 
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If  now  /a  this  suggests  that  the  sources  remain  zero.  This 

must  certainly  be  true  for  physically  meaningful  solutions, 3  » 
for  the  magnetic  field;:  if  no  magnetic  field  is  present  originally, 
none  can  be  generated  spontaneously  by  the  fluid  motion. 

To  substantiate  analytically  t his  p hysical  interpretation 
of  (5)  one  must  inquire  into  the  nature  of  the  solutions  of  the 
integral  equation  to  which  (5)  leads.  Prom  (4-14)  we  have 

fctr**)  =  -v"1|df'G(T,t  |  t\Y)7T(  tW  )-■ 

+v"1jd3r/G(  r,t\r',o)  bir',0) 

/  /  3  / 

in  which  d  P  =  dt  d  T  ,  G  is  the  appropriate  Green’s  function, 

T«J 

is  by  (4-11)  a  solution  of 


(6) 


where 


TT( f\t)  =  A(r,t)  *  d r.t  \  T,t  )7T(  r#t  > 

J 

'  A(T.t)  =  Jd3r'/<(r,t|  rlo)  b(rio) 


(V) 


K(  T,t|  rlt'i  =  -v“V.7cKrrt|r'.t'>  . 


-i 


/  / 


( 8^ 


is  to  be  regarded  as  a  known  function.  In  what  follows 
|l  =  T/(  T-  pt  )  etc.,  and  ^  will  denote  the  combination  of 
variables  ('f\t). 

One  approach  to  solving  (7)  is  ■  the  method  of  iterations 
(Page„  sec.  9.2;  Co urant »  Hilbert,  Ch.  3.  sec.  6).  Choose  any 
finite  continuous  vector  function  as  the  zero  order  approximation 
for  7 r  ;•  let  it  be  called  T  •  Then  the  first-order  approximation 
for  ~fj~  is 


=  Acf  >  +jdj/K(fipjrty)t 
Tir  a<s>  K(^|f')  A</> 

+|ds>/  I?) 


(9) 


(10) 


i  L. 

:  f 


*  !•:• 

I  «  w 


..  . 
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where 


K2(y/s')  • 


If  this  iteration  process  is  continued  indefinitely  and  we  define 
the  iterated  kernels 

n^2  (1. 

K.<flp  -Kifls'y. 

there  obtai  ns  for  the  n—  approximation 

Tn<P>  =  A<^>  *Z  A(/) 

*S,J  (i: 

It  is  shown  in  treatises  on  integral  equation  (Whittaker  and 

Watson,  1920,  p.  222)  that  regardless  of  the  choice  of  a  trial 
solutlc  T°  ,  the  last  t erm  in  (13)  approaches  zero  and  the 
series  converges  uniformly.  Hence,  the  solution  of  (7)  may  be 
had,  in  principle,  to  any  desired  degree  of  accuracy.  If  we^ 
moreover,  assume  for  simplicity  that  the  series  of  kernels 

L(f  |/)  =  X  K^|f;)  d' 

t * 

also  converges  uniformly,  the  expression  (j 3)  can  be  put  in  a 
concise  forms 

Tf(j)  =  A<$)  +  L(^  |^)  A(/)  •  a* 

Li  is  called  the  reciprocal  kernel. 

A  glance  at  (8)  shows  that  if  b(  T»°)  =  0  then 

=  0,  and  from  (6)  there  follows  b  (  T,t)  =  0.  The 
physical  requirement  that remain  zero  if  initially  zero  is 
met  by  our  solutions,  and  this  is  true  of  divlB  as  well: 


if  7*  8  (  T»0)  =  0  then  V*  -B  (  T,t)  =  0. 


(16) 
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B.  (Condition  (3),  is  not  in  general  satisfied  by  solutions  B 

of  equation  (l)j:  at  least  it  is  not  necessary  to  impose  (3)  on 
the  quantity B  appearing  in  the  partial  Integral  (3-2).  As  a 
result  the  non-homogeneous  term  in  (4)  must  be  augmented  by  an 
additional  termfrom  (2),  giving 


(17) 


|g(  v-te  )  +  V*  V(  V*  £  )  +v*(  V*  V  )£>  =  v72( 7«B). 


Repeating  the  preceding  steps  for  this  equation  shows  that  if 
<  v-ni-o,  then s  nevertheless,  V- (  \/ •  V)B  develops  into 
a  source  of  ,  even  if  Q  V  *  P7  )£Gt_Q=0« 

Solutions  of  the  substantial  diffusion  equation 
are  divergence-free  if  and  only  if  ^7*  ( V*  7YB  = 

For  the  partial  integral  under  the  circumstances  of 
sections  3d  and  3e ,  however,  it  is  not  essential  that 
d  =  0,  but  rather  that 

V°*  £  =  0,  (18) 

by  virtue  of  (3-13) .  Consider  the  substantial  diffusion 
equation  (1)  in  the  form 

A 

af=v-V2B.  as) 

/V  A  _  0  — 1  0 

if  B  is  expressed  in  Lagrangian  terms,  B[r(rc  ',t),tj  =B(r  ,t), 
then  (19)  may  be  written 

|f  =  vy2  £  •  (10l> 

d/8t  and  ^7°  non  commute  and 

^(V°-B)  =  V  V°.  \72B  •  (20) 

To  analyze  this  equation  the  mixed  derivatives  on  the  right 
must  be  replaced  by  derivatives w  ith  respect  to  the  Lagrangian 

O  2  A 

coordinates  T  •  ^  B  may  b®  expressed  as 


36 


a  Bi 

=?*j 


a 

.  < 

A 

3B. 

sj 

3*5’ 

P  ^ 

3dB. 

M- 

-.2  o  ^ 

3  xk  “i 

sxj 

3x .  3x,° 
j  k 

3xj9x3  3xk 

jfS 

3x°  a2^ 

-.2  o 

3  xk 
-H-  K 

3^i 

ax3 

ax  j  3x 

0-.  0 
i^k 

dx]dk3 

34 

(21) 


W’e  consider  first  the  linear  velocity  shear  of  section 
3c.  The  Lagrangian  trajectories  are 


x-,  =  x°  ,  x0  =  x°  +  (<2.x°  +.  b)tr,f,  x,  =  x°. 

-L  X  d  d  o  o  o 


(22) 


after  (3-18)  and  (3-26)  j  thus. 


V2r°  =  v°2f  =  o 


(23) 


which  reduces  (21)  to 


O 

a^B, 


axjaxj 


=  (V^*7 


2^ 

O  a  b4 


'll 


3x.  3x, 

k  t 


(24) 


For  7°  «  v  B  we  have 


33S 


1  =(Vx°-C7x“) 


_3^ 
3  B4 


3x.  3x  .3x  . 
i  J  J 


roTo,  o 

3xk3x«  3xi 


(25) 


32B 

*— —o  4  ^4  -V4)  ■ 

Kd4  K  1 


The  following  relations  are  easily  obtained  from  (22) 
Vxx  =^i»  ^x2  =  ^*3  =^3  * 


(26) 
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7*°  .  Vxl  =  7*°  .  =  1 

7x°  .  7*°2 , 7^  •  w3  «  o 

Vx°  .  V4  =  1  *^2t2 
Vx°  •  7x3  =  -  at. 


(27) 


The  last  term  in  (25)  vanishes,  leaving 


v°  -VZB  ^  V°  ( V°-  B  )  -  2  at  -3^-0  (  V  B  > 

2  3 

t  V°-B  >, 


+  a  t 


2,2  S2 


3'(x°)2 


and  (20)  becomes 


=  V°2(V°-B)  -  2a.t  (V°-B> 

3x23x3 


(28) 


+  it!  -2-5-8  (  V°-8>- 


S(x2) 


(29) 


t  = 


Note  that  in  general  d(  (7  •  B)  /dt  p}  however,  when 

t 

0  the  equation  is  momentarily 


^(7°-B)  =  V  V°2(7°vB). 

As  before,  we  see  that  if  V°*B  =  0  initially,  there  is  no 
tendency  for  sources  to  develop,  and  we  surmise  that'  V 
remains  thereafter  zero.  This  inference  is  given  added  w eight 
by  considering  what  happens  when  t  is  small  but  not  zero. 

The  last  term  on  the  right  of  (29)  can  be  ignored  compared 
with  the  first,  and  for  t->0, 
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-v"1  &  (V°'B)  =  2«t  -4-0(  V°-B). 


3x“3x” 


For  the  sake  of  brevity  we  define  the  operators 
,S  _  3S  _  aS 

3Xg3x°  ’ 

and  let  j)  designate  the  combination  of  variables  (  T°»t)  while 

p(f)  =  7°*  B  • 

Equation  (3)  in  this  symbolism  becomes 

V°  D  -v”1^  =  2  0,t32D.  (31.) 

Comparing  this  with  (4-1  )  and  (4-14)  we  see  that 

D(j>)  =  Jdf/G(  J  |^)2at9/2D('/)  +  v*’1Jd3T°/(GD/)t  =0  •  (52) 

Operating  on  this  with  with  32  gives  an  integral  equation  for 
d2D,  viz.. 

Jf  2 

d/2at32Q(^)S  D(f  ),  (33> 


where 


(j>)  =  =o  • 


It  has  a  solution  of  the  form 


3?D(j>).  =  Cpif)  •*Jdj>/M(j>| />^(/)» 

after  (IV-7)  and  (IV-15),  M  being  the  kernel  reciprocal  to 
2&t92Cr.  As  before,  if 


Dt=0  =  (7°‘£(r°,o)  =  o 


then  ^(f)  =0  by  (34)  ,  and  consequently  32D  =  D  =  0  by 
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(35)  and  (32). 

Interpreting  this  result  is  not  quite  so  simple  as  before. 

D  is  expressed  in  Lagrangian  terms,  so  the  partial  time  derivative 
in  (31)  is  equivalent  to  a  substantial  derivative;  accordingly, 
what  we  have  shown  is  that  if  D('/'O,0)  =  0  then,  subsequently, 

D(  T°t)  =  0  at  a  particle  moving  with  the  fluid.  But,  since  this 
must  be  true  for  all  particles,  we  conclude  that 

7  °  *  £  ( T°t )  =  0  everywhere  in  the  fluid  and  at  all 
times  provided  only  7°*3(  T^O)  =0  (37) 

Equation  (3-14)  follows  at  once. 

As  regards  the  approximate  solutions  of  section  3d, 
the  last  terms  of  (21)  and  (25)  are  no  longer  zero.  Nevertheless, 
since  V  is  slowly  varying,  the  curvature  of  the  trajectories, 
of  which  the  second  derivatives  of  the  Lagrangian  coordinates 
in  (21)  and  (25)  are  a  measure,  is  very  small,  and  the  last 
terms  In  these  equations  are  negligible  compared  with  the 
first  terms  and  may  be  dropped,  leading  once  again  to  (29). 

G.  It  remains  to  treat  the  divergence  criteria  for  the 

solutions  of  section  3c.  The  essential  steps  are  exactly 
the  same  as  in  A  of  this  appendix.  Taking  the  divergence  of 
(4-18b)  shows  that 


V2(V-F)  +  k2  V-F  =  v"1'/.  V(  V*F)  (.38) 

provided 

7°  (  P*  V  )  V  =  0  .  (39) 

v  v  y 

Thus  we  let  f  represent  either  p  or  V  mF  and  /  either  P  or  7*p'. 
Then  from  (4-24)  and  (4-25)  , 

f  ~V~1  (40> 


(41) 

(v*7>f  =  (v-V)f  d5r/(v-7G>(v*  7)/^ 


and 
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which.  have  precisely  the  same  forms  as  (6)  and  (7),  respectively. 

We  may  state  at  once  that  if  $  =  0  then  both  (V*  V )  £  and.  £  vanish. 
*  Once  a  particular  solution  of  the  homogeneous  equation  is  selected, 

p1  is  given  by  (40),  and  this  fixes  the  form  of  V  *  F  t.  it  is 
that  for  which  the  "inhomogeneous"'  term  is  V*  P  .  And 

v 

if  V'F  -  0  then  V*  r  =0.  (42) 

This  leads  to  equation  (4-30)  . 

D.  From  the  physical  and  mathematical  analyses  of  the 

'"inhomogeneous"'  differential  equations  studied  in  this  appendix, 
we  hazard  a  generalization.  Consider  an  equation  of  the  form 

V2A  = 

where  A(Tj'b)  is  an  arbitrary  vectory  function  and  denotes 
an  operator  linear  in  the  space  derivatives.  If  /\ (  T% 0)  =  0 

we  assert  that  since  £A  (7^,0)  =  0,  a  result  corresponding  to 

■»  A 

(15),  there  are  initially  no  sources  for  A  ,  and  moreover, 

because  of  this  there  is  no  tendency  to  generate  sources. 

Whenever  the  latent  "sources,"’  -if  A,  depend  on  A  in  such  a 
way  that  they  are  zero  when  A  everywhere  zero,  they  are 
given  no  chance  to  develop,  and  the  quantity  A  *  once  zero 
remains  zero  thereafter. 


References 


-41- 


Brand,.  L.,  1947.  Vector  and  Tensor  Analysis  (New  York:  Wiley), 

p.  263. 

Courant,  R. ,  and  Hilbert,  D.,  1953.  Me thods  of  Mathematical 

Physics „  Vol.  I.  (New  York;  Interscience). 

Elsasser,  W.,  1954.  Phys.  Rev.,  95,  1. 

Inglis,  D. ,  1955.  Rev.  Mod.  Phys.,  27,  212. 

Lundquist,  S.,  1951,.  Phys.  Rev.,  83,  .307. 

1952.  Arkiv  Fysik,  _5,  297. 

Morse,  p.,  and  Feshbach,  H.  s  1953.  Methods  of  Theoretical 

Physics .  Vols,  I  and  II.  (New  York:  McGraw-Hill). 

Page,  C . ,  1955.  Physical  Mathematics.  (Princeton:  Van  Nostrand)„ 

Skabelund,  D. ,  1955.  Technical  Report  17  of  this  series. 

University  of  Utah,  Salt  Lake  City, 

Sommerfeld,  A.,  1949.  Partial  Differential  Equations  in  Phys ics . 

(New  York:  Academic  Press) . 

Whittaker,  E. ,  and  Watson,  J. ,  1920.  Modern  Analysis ,  3d.  Ed. 

(Cambridge  Univ.  Press). 


-42 


DISTRIBUTION  LIST 


Addressee  No «  of  Copie 

U.S.  Navy  Hydrographic  Office,  Washington  25,  D.C.  1 

Director,  Woods  Hole  Oceanographic  Institution, 

Woods  Hole,  Massachusetts  1 

Director,  Scripps  Institution  of  Oceanography. 

La  Jolla,  California  1 

Director,  Lamont  Geological  Observatory 

Torrey  Cliff,  Palisades,  New  York  1 

Director,  U.S.  Geological  Survey,  Department  of  Interior, 

Washington  25,  D.C.  1 

U.S.  Geological  Survey,  Geophysics  Branch, 

Washington  25,  D.C.  1 

Director,  Institute  of  Geophysics,  University  of  California 
Los  Angeles,  California  1 

Princeton  University,  Department  of  Geology, 

Princeton,  New  Jersey  1 

Carnegie  Institution  of  Washington,  Department  of 

Terrestrial  Magnetism,  Washington  25,  D.C.  .  1 

Dr,  G.P.  Wollara,  University  of  Wisconsin,  Madison, 

Wisconsin  1 

Geophysics  Branch,  Code  416,  Office  of  Naval  Research, 

Washington  25,  D.C.  2 

Director,  Naval  Research  Laboratory,  Attention: 

Technical  Information  Officer,  Washington  25,  D.C.  6 

Off icer-ln-Charge ,  Office  of  Naval  Research,  London  Branch 
Off..  Navy  No.  100,  Pleat  Post  Office,  New  York,  New  York  2 

Office  of  Naval  Research  Branch  Office, 

346  Broadway,  New  York  13,  New  York  1 

Office  of  Naval  Research  Branch  Office,  Tenth  Floor, 

The  John  Crerar  Library  Building,  86  East  Randolph  St. 

Chicago,  Illinois  1 

Office  of  Naval  Research  Branch  Office,  1030  East  Green  S't. 
Pasadena  1,  California  1 

Office  of  Naval  Research  Branch  Office,  1000  Geary  Street,  ■ 

-  SanFrancisco,  California  1 

Office  of  Technical  Services,  Department  of  Commerce, 

Washington  25,  D.C.  1 

Armed  Services  Technical  Information  Center,  Documents 

Service  Center,  Knott  Building,  Dayton  2,  Ohio  5 

Assistant  Secretary  of  Defense  for  Research  and  Development, 
Attention:  Committee  on  Geophysics  and  Geogrg?  hy. 

Pentagon  Building,  Washington  25,  D.C.  1 


